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q_J Abstract: We continue to study the holographic QCD (hQCD) model, proposed in a 

previous paper, in an Einstein-Maxwell-Dilaton (EMD) system. In this paper we discuss 
some aspects of quark gluon plasma (QGP) in the hQCD model, such as drag force, jet 
quenching parameter and screening length. The results turn out to be consistent with those 
as expected in QCD qualitatively. By calculating free energy of the background black hole 
solution, we find that there exists a Hawking-Page phase transition between small black 
hole and big black hole when chemical potential /i is less than a critical one /x c , and the 
phase transition is absent when chemical potential is beyond the critical one. 
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1 Introduction 

The analysis of various dynamical quantities from the experimental data obtained at Rel- 
ativistic Heavy Ion Collider (RHIC) in Brookhaven National Laboratory leads to strong 
evidence of appearance of strong coupled quark gluon plasma (QGP), a deconfined phase 
of QCD at high temperature and high number density [1], [2], [3]. In the Au+Au collision 
with the maximum center-of-mass energy around 200 GeV, several phenomenological fea- 
tures, e.g, a very small value of shear viscosity, quenching of high energy partons with large 
transverse momentum, and elliptic flow etc, indicate that the dynamics of thermal medium 
produced after collision is dominated by non-perturbative effects [4]. Being formulated in 
Euclidean time, lattice QCD seems to be the best candidate to explain the phenomena in 
thermal equilibrium, while the perturbative QCD works only in the weak coupling region. 
Both the perturbative QCD and lattice QCD are failed to compute some dynamical quan- 
tities like transport coefficients, drag force, and jet quenching parameter etc. Due to the 
feature of strong/weak coupling duality, the AdS/CFT correspondence [5] [6] provides a 
powerful tool to stress those issues. 

In the AdS/CFT correspondence, the well-known example of strong/weak duality is 
AdSs/CFT4 and in this case, the four dimensional conformal field theory is J\f = 4 SU(N C ) 
super Yang-Mills (SYM), while the bulk theory is the type IIB supergravity in AdSs. In 
the AdSs/CFT4 correspondence, the ratio of shear viscosity over entropy density for the 
dual field theory is found be small (r//s = l/4ir) for large number of colors (Nq) and large 
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't Hooft coupling (A = gYM 2 Nc) [7] [8] [9]. It turns out that this value is consistent with 
the one from RHIC data with experimental error [10]. In addition, it is found that the 
value 1 /4tt is universal for various deformations of the gauge theories with gravity duals. 
Another example of the universality is discussed in [11], where it is found that for some 
gauge theories (not necessarily conformal) with well defined gravity duals where if some 
conditions on the bulk stress energy tensor are satisfied, the electrical conductivity at finite 
chemical potential (//) and temperature (T), the thermal conductivity (kt) and the ratio 
of thermal conductivity to viscosity fi 2 ) are independent of any specific model. 

Certainly a realistic holographic model dual to the strong coupled QCD at finite tem- 
perature and finite density is a good starting point to study the RHIC physics. Unfor- 
tunately such a model is still absent. However, due to the universality mentioned above, 
one is expecting to understand various features of QGP at RHIC with some deformed 
AdS^/CFTi, such as drag force, jet-quenching parameter, and screening-length etc. In- 
deed, over the past years a lot of works have been done along this line. 

By employing the so-called potential reconstruction approach [12-15] to build up holo- 
graphic QCD model [16] from bottom up point of view, they proposed a hQCD model in 
an Einstein-Maxwell-dilaton theory in [15]. In [14, 15], some properties of the hQCD 
model are studied, such as equation of state, Wilson line operators [12] and confine- 
ment/deconfinement phase transition and associated phase diagram [15], with good agree- 
ment with lattice QCD simulation. It shows that this model can capture some most im- 
portant characteristics of real QCD. In this paper, we continue to study the hQCD model 
by investigating some aspects of QGP phase such as drag force, jet quenching parameter as 
well as screening length. The results show that these quantities are consistent with QGP 
properties. In this sense, we further confirm that this model may give us the hints on the 
efficient way to deform the pure AdSs geometry and to realize the holographic description 
of low energy QCD. 

When high-energy partons perform a dragged motion as they pass through the QGP 
medium, their energy loss can be encoded by drag force. In the AdS/CFT correspondence, 
an external heavy quark and its gluonic neighborhood are mapped into the endpoint of 
fundamental string attached to the AdS boundary and the string itself in the AdS bulk 
geometry, respectively. This external quark, with mass proportional to the length of the 
string, loses its energy as the string trailing back imposes a drag force on it. Within the 
framework of gauge/gravity duality, the drag force experienced by an external heavy quark 
moving with a constant velocity in N = 4 super Yang-Mills plasma at finite temperature 
is computed [17], [18], [19], [20], [21], [22], [23]. There are also further generalizations of 
drag force computation for the charged N = 4 SYM [24] and with the backreaction effect 
due to the static heavy quark cloud distribution [25] . 

Due to the medium, the suppression of heavy quark with high transverse momentum 
leads to energy loss which is so called jet quenching phenomenon [26], [27], [28]. The 
transport coefficient q [29], characterizing such phenomena, is defined perturbatively as 
the ratio of square of the mean transverse momentum over the mean free path [30] . With 
the framework of eikonal approximation [31], [32], the parameter q can be calculated from 
light-like Wilson loop in adjoint representation [33], [34], [35], [36], [37], [38], [39]. Ref. [40] 
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and [41] calculate the jet quenching parameter in the presence of chemical potential. Using 
AdS/CFT, the quark-antiquark pair is mapped to the two endpoints of a fundamental 
string, both endpoints attach on the AdS boundary. The string configuration shows it hangs 
down to the bulk along radial direction and turns back to the boundary again. The light-like 
Wilson loop in the fundamental representation of boundary gauge theory and its thermal 
expectation value correspond to the trajectory of two endpoints of the string and exp(<S) 
respectively, where S is the string Nambu-Goto action. In addition, we further extend 
our analysis by studying the time-like Wilson loop and relate it with another important 
parameter, the screening length of a quark-antiquark pair. It is defined as the maximum 
length between a moving qq pair, beyond which they break off with no binding energy 
and thus become screened in the QGP medium. Screening length depends on the velocity 
and the orientation of the quark-antiquark pair with respect to the medium. Treating 
J\f = 4 SYM at finite temperature as a boundary theory, the binding energy between quark 
and antiquark pair moving in the QGP and the screening length are calculated in [42]. 
There are also other important generalizations to calculate these quantities by introducing 
a boost on the original static background [43], [44], [45]. In this paper we will calculate 
the screening length in the static frame of qq pair for our hQCD model. 

In [15] we calculated heavy quark potential between a quark-antiquark pair in our 
model, and found that there is a confinement/deconfinement phase transition in the hQCD 
model, and that there is a critical point in the T — fj, phase diagram. In this paper we 
further confirm this phase transition in the small (i region by computing free energy of 
the background black hole solution by using the method in [46]. In the AdS/CFT corre- 
spondence, the Hawking-Page phase transition [47] between AdS black hole and thermal 
gas in AdS space is identified with the confinement /deconfiement transition in gauge the- 
ory [48] [49] [50] [51]. In our case, the phase transition happens between small black hole and 
big black hole, which will be clear shortly. In addition, we will argue that there is no phase 
transition between black hole solution and thermal gas solution in our model. 

The organization of the paper is as follows. In section 2 we briefly review the poten- 
tial reconstruction approach to the Einstein-Maxwell-Dilaton system by generalizing the 
discussion in [15] to the case with a coupling between dilaton field and Maxwell field. In 
section 3, we discuss the generic black hole solutions with asymptotical AdS boundary, and 
in particular present an analytic black hole solution. In addition, in this section we also 
briefly review the black hole solution for the hQCD model studied in [15]. In section 4, 
we calculate the drag force in this hQCD model. Jet quenching parameter and screening 
length are discussed in section 5 and 6, respectively. In section 7, we study the free en- 
ergy of the background black hole solution and discuss the phase transition between small 
black hole and big black hole in small chemical potential region. Section 8 is devoted to 
conclusions and discussions. 

2 Einstein-Maxwell-Dilaton system 

In this section, we use the potential reconstruction approach [12, 15] to study a 5D Einstein- 
Maxwell-Dilaton (EMD) system. In [15], the authors did not consider the coupling between 
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gauge field and dilaton field in Einstein frame. Here we take the coupling into consideration 
in a more general version 

S 5D = J d^v^Fe- 2 * (r s + 43 M 0#V - V S (4>) - ^-e^F^F^ , (2.1) 

where the action (2.1) is written in string frame, F^ = d^A u — d u A^ is the Maxwell 
field, Z{<p) is an arbitrary function of dilaton field <f> and V s ((f>) is the dilaton potential. In 
Einstein frame we can rewrite the action as [12] 

S,d = J tfxyf^ (r - ~ Ve(<P) ~ ^ V") , (2.2) 

-4<f> 

where Vs = Ve& 3 • The metrics in these two frames are connected by the scaling trans- 
formation 

9% = e%&. (2.3) 
The Einstein equations from the action (2.2) read 

E,u + \gfL + V E {^ ~ \d^dA - ^ (f^fS - \g* v F u F u ^ = 0(2.4) 



where E^ v = R ilv —\Rg iiV is Einstein tensor. We here consider the ansatz A = Ao(z)dt, 4> = 
(f)(z) for matter fields and 

ds 2 s = (-f(z)dt 2 + y^y + dx'dxA , (2.5) 

for the metric in string frame, where i = 1,2,3, I is the radius of AdSs space, and A s is 
the warped factor, a function of coordinate z. The metric in the string frame will be used 
to calculate the loop operator below. In Einstein frame the metric reads 

£ 2 e 2Ae ( dz 2 
ds\ = 5— -f(z)dt 2 + — — + dx % dx l 



z 2 V W /(*) 



e 2 p 2A s -f / j 2 

-f(z)dt 2 + ^+ dx { dx l ) , (2.6) 



z 2 v w m 

with A e = A s — 20/3. In the metric (2.6), the (t,t),(z,z) and (xi,Xi) components of 
Einstein equations are respectively 

b'(z)f(z) b\z) 2 4 2 A o i(z) 2 Z(0) V E (<P)b(z) 2 n 

2 _ %'{z)f{z) _ 9b'(z) 2 _ 34/(*) 2 ^) _ W E {4>)b{z) = n 
9U 8b(z)f(z) 4b(z) 2 8g 2 g b(z)f(z) Af(z) 

, 3b'(z)f(z) 4 , 2 3/(z)6"Qs) 3/(z)fr'(^) 2 Mg)^) T , 

7 + b(z) + 3 /( ^ (z) + ~ &£) Mz) 2 2gjbW~ + Eim ) = ' 

(2.7) 
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where b(z) = £ 2 e 2Ae /z 2 , and Aq(z) is electrical potential of Maxwell field. From those 
three equations one can obtain following two equations which do not contain the dilaton 
potential Ve{<P), 



A'*{z) + A' s {z) 



40V) 2 
3 z 



f"(z) + f'(z)[3A' s (z)-2 ( f>'(z) 



A' s (z 
3 



2 _ 2f(z) 44f(z) 
3 

40(z) 



0. 



3 3z 
z 2 Z(<j>)e^- 2A °MA' (z) 2 



9lL 2 



0. 



(2.8) 
(2.9) 



Eq.(2.8) is our starting point to find exact solutions of the system. Note that Eq.(2.8) 
in the EMD system is the same as the one in the Einstein-dilaton system considered in 
[13] [14] and the last term in Eq.(2.9) is an additional contribution from electrical field. In 
addition, the equation of motion (EOM) of the dilaton field is given by 



:d z 



0: 



^W- T ^ w) + ^)^ W a =0| (21Q) 



And the EOM of the Maxwell field is given as 



V- 



(2.11) 



From equations of motion, once A s (z) is given, we can obtain a general solution to the 
system, which takes the following form 



z e 2Mx) (1 lo y 2 e- 2A ^A' s (y) 2 dy + <h) 3A s (z) 



Aq{z) = A Q0 + A 



01 



ye 



2 -^-A s {y) 



V E (z) 



a .3 e 2</,(x)-3A s (x) 



( ^ do 



dy 



(2.12) 
(2.13) 



Z{<t>{y)) 



dy 



9V 2 



+ fi 



\ 



dx + f , (2.14) 



,- 2A ^^z 2 f{z)j e- 2A ^ +iJt ^Z{<P(z))z 2 A' {z) 



Ag 2 £ 2 f(z) 

2 (3 + 3z 2 A' s (z) 2 + Az4>'(z) + z 2 4>'{z) 2 - 2zA' s (z) (3 + 2z<f/(z))) 



f(z)(-3 + 3zA' s (z)-2z^(z)) 
2zf(z) 



(2.15) 



where (j>o, Aqq, Aqi, fo, fi are all integration constants and can be determined by suitable 
UV and IR boundary conditions. When Z(cp) = 1, the general solution reduces to the one 
given in [15]. Thus we have given a generic formulism to generate a set of exact solutions 
of the EMD system with a given A s (z). 
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3 General asymptotical AdS black hole solutions 



Since here we are only interested in the black hole solutions with asymptotic AdS boundary, 
we impose the boundary condition /(0) = 1 at the AdS boundary 2 = 0, and require 
(p(z), f(z), Aq(z) to be regular at black hole horizon Zh and AdS boundary z = 0. There 
is an additional condition Ao(zh) = 0, which corresponds to the physical requirement that 
A^A^ = g tt A A must be finite at z = Zh- 

We can parameterize the function f(z) in Eq.(2.12) as 



12 



01 



(/o* 9(r)dr J: 9 -^)dx 



Jo g(x)dx 
Io h g{x)dx ' 



where /o 



A 2 

1 f , — oi 
L > J 1 ~ 4gjF 



/o h 5 (a;)dx 



and 



(3.1) 



(3.2) 



We expand the gauge field near the AdS boundary to relate the two integration constants 
to chemical potential and charge density, respectively, 



with 



A {z) 



A 00 
An 



e -3- A a (y) 

Aqo + Aqi NN z 2 -f 



z(<P(y)) 



fj- 



/' 



An, f Z h 9{y)\ J.. 

zwy)) ay h z{mr v 



(3.3) 



(3.4) 
(3.5) 



The temperature of the black hole can be determined through the function f(z) in 
(3.1) as 



1 

4-7T 



!/'(*)! 



z = z h 



A 2 

^01 



g(z h ) ^ g{r)dr ^dy - g(z h 



fo h 9{x)dx 



(3.6) 



Following the standard Bekenstein-Hawking entropy formula, from the geometry given in 
Eq.(2.6), we obtain the black hole entropy density S as 



S 



A n 




2/. 



(3.7) 



where V3 is the volume of the black hole spatial directions spanned by coordinates X{ in 
(2.6). 
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3.1 An analytical black hole solution 

In this subsection, we list an analytical solution of the Einstein-Maxwell-Dilaton system 
by using Eq. (2. 12-2. 15) with Z((p) = 1. We impose the constrain /(0) = 1, and require 
4>(z), f(z) to be regular at z = 0, and z^. We give the solution in Einstein frame 

/>2 2A e / l 2 \ 

da| = iiL_ ( -/(z)^ 2 + — + dzWj , (3.8) 

with 



^e(z) = log 



z 



£ sinh(^; 



/(z) = 1 - — ii(3sinh 4 (-) + 2sinh b (-)) + -V^ 2 sinh 4 - 

3 Zq Zq 8 \ Z 

- h 

A (z) = ^V 12 stall 2 (3.9) 

where zo is an integration constant and Vu, V\i are two constants from the dilaton potential 

12 + 9 sinh 2 ( + 16Vh sinh 6 ( t) V? 2 sinh 6 ( ^■ 



= ^ ^ + " 8^ - (3 - 10) 

The two integration constants V\\ and V\i then can be expressed in terms of horizon z^ 
and chemical potential [i as 



3cosh 4 ( | ^ 

V U = — 



\2zo) I ~~ 4^ 



9 




M^ocosh 



32 (2sinh 2 (^)+3 

2 f 



We can obtain the temperature of the black hole by using Eq.(3.6) based on the above 
formulas. In Figure [1], we show the temperature as a function of horizon radius z^ in 
cases of three different chemical potentials /x. In this plot we take parameters I = l,zo — 
l,g g = 1. We see from Figure [1] that the temperature with respect to horizon Zh is 
monotonic for a fixed chemical potential. A vanishing temperature means that the black 
hole is extremal with a smallest horizon radius. The smallest horizon radius increases as 
chemical potential becomes large. 

3.2 The hQCD model 

Based on the general solutions, in Ref. [15], a holographic QCD model is proposed to realize 
the confinement /deconfinement phase transition of QCD. Since the aim of this paper is to 
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Figure 1. The temperature as a function of horizon radius z% for the analytical black hole solution 
with parameters t = 1 , zq = 1 , g g = 1 . 



further study the hQCD model, we here briefly review the main features of the model. In 
our hQCD model, the warped factor A s (z) takes the form 

A s (z) = k 2 z 2 , (3.12) 

where A: is a constant. To take this factor has various phenomenological motivations in 
order to build a successful holographic QCD model, for details see [15]. In this paper, we 
set k = 0.3GeV, following Ref. [15]. With this factor, we have the dilaton field <p as 

4(z) = \k 2 z\l + H{z)), (3.13) 
where we have set the integration constant <f>Q = 0, and H (z) is given by 

H{z) = 2 F 2 (l,l;2,^;2£; 2 /) . (3.14) 
The characteristic function of the black hole background takes the form 

Jo 9i x ) (fo h 9{r)dr g(y)*dy) dx 



m= 1 + 



/ \ 2 r z 

1 / n 



^ 2 \f Zh g(y)Uy) 9{x)dx 
Jo 9{x)dx 
Jq H g{x)dx' 



(3.15) 



where 



g(x)=x 3 el k2x2 ( 1+H ^- 3k2x \ (3.16) 



One can clearly see that the second term in (3.15) comes from the contribution of electric 
field. If one turns off the electric field, one can reproduce the black hole solution in Einstein- 
dilaton system [13]. In addition, the electric field At(z) is given by 

A t { z ) + f xe^^~ 1+H ^dx. (3.17) 

J Zh g(y)idyJo 
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Figure 2. Plot (a): The black hole temperature as a function of horizon with different chemical 
potentials. When \i > \x c the temperature monotonically decreases to zero with increase of z^; when 
< fi < fx c , the temperature decreases to a minimum at z m and grows up to a maximum at zm and 
then decreases to zero monotonically. When /i = [i Cl one has z m = zm- Plot (b): The temperature 
of the black hole with = O.lGeV. The three black hole solutions with horizon z' p , z p and z' p ' have 
the same temperature. The black hole with z m < z p < zm is thermodynamically unstable. Here 
we take g g £ = 1, k = 0.3GeV. In this hQCD model, we always fix k = 0.3GeV and accordingly the 
critical chemical potential is \x c = 0.34GcV, which corresponds to the case z m = zm- 

The temperature behavior of the black hole solution is discussed in some details in [15], 
with respect to horizon radius and chemical potential. To be complete and for later use, we 
here briefly repeat some main salient features. In Figure [2] we plot the temperature with 
respect to horizon radius with different chemical potentials. One can see clearly that 
the temperature behavior crucially depends on the value of chemical potential: there is a 
critical chemical potential /x c , beyond which the black hole is always thermodynamically 
stable, while when the chemical potential is less than the critical one, there is a region 
of horizon radius, where the black hole is thermodynamically unstable with negative heat 
capacity. To be more clear we plot in Figure [2](b) the temperature versus the horizon Zh 
in the case \i = O.lGeV < [i c as an example. One can see from the figure that the black 
hole is thermodynamically unstable in the region z m < Zh < zm , where z m and zm are the 
black hole horizons corresponding to the minimal and maximal temperatures, respectively. 
In this region, the heat capacity of the black hole is negative. The black hole solutions in 
the regions z^ < z m and z^ > zm are thermodynamically stable. When fj, > fj, c , z m and zm 
are degenerated to one point. In contrast to the case in Figure [1], there are local minimal 
and maximal values of temperature in small fi cases. This is crucial to realize the critical 
point in the T — \i phase diagram of the hQCD model [15]. 

4 Drag force 

In this section, mainly following [18] we compute the drag force experienced by an external 
probe quark traversing in the QGP in the present hQCD model. The drag force is a function 
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of temperature and chemical potential. The boundary gauge theory we are considering is on 
A4 4 described by the boundary coordinates t . X - X - X . The dynamics of a fundamental 
string is completely specified by the Nambu-Goto action in the black hole background 
within the string frame (2.5). On this background the world sheet action reads 

1 f dX^ dX^ 1 

S = -^J drdaV-det 9a „ = ——G^ (4.1) 

where g a p is the induced metric on the world sheet and G^ v is the background metric. The 
equation of motion derived from (4.1) is given by 

A a P% = 0, P« = - 1 ±- i G ilv d<»X\ (4.2) 

where A a is the covariant derivative with respect to g a p and P° is the world sheet current 
of space time energy-momentum of the test string. We consider the motion of the string 
along x . In the gauge, r = t and a = z, the string dynamics can be completely specified 
by the function x (t, z). In this case, the Lagrangian reads 



11 f(z)(d z x^ (dtx^l 

+ H Wiz)' ( } 



where H is defined as 



z 2 



To capture the dragged motion of the quark in the boundary theory we assume the following 
ansatz in the bulk [18] 

x l {t,z) = vt + i(z). (4.5) 

Here we have assumed only the late time behavior of the string motion. With this ansatz 
the Lagrangian reduces to 



1 / 1_ + f{z){d z g Z )f _ 



2na' \j H H Hf(z 

The momentum which conjugates to £(z) reads 



n 9C g / / Hf 

For the sake of consistency it is important to invert the equation (4.7) and write it in the 
following way 



n f 2 (/-^ 2 ) , x 

(4.8) 



Here the positive sign is taken due to the trailing nature of the string profile [18]. To obtain 
the string profile we have to solve the differential equation (4.8). To have a real £(z), we 
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further impose the constraints 



mi 



lie 



1 



so that one has £'\ z=Zv = v 2 / f 2 , keeping finite. The profile of the string is defined in the 
region with z < z v , That is, there is a maximal value z v < Zh for the string profile. 

The constraints are very useful to figure out the final form of the drag force. Before 
to compute the drag force, we here mention the relation between the drag force in the 
boundary field theory and the dissipation of momentum flowing down the string, in light 
of AdS/CFT correspondence. In the boundary theory the presence of the thermal medium 
results into dissipation of energy and momentum of external quark until it reaches thermal 
equilibrium with the medium. In the bulk theory the momentum is flowing down the string 
from the boundary to the bulk and the change of momentum at a given spatial point on the 
world sheet for a given time interval can be calculated. The identifications of the endpoint 
of the string attached to the boundary with the quark and of the string in the bulk with 
the thermal medium around the quark suggest that the drag force can be realized in terms 
of the force imparted by the string on its boundary endpoint. To calculate the change of 
string momentum due to its motion along x\ direction, we consider a closed curve on the 
world sheet and study how the momentum is conserved around this curve [55] . According 




Figure 3. This plot shows a closed path in an anti-clockwise direction on a world sheet bounded 
by coordinates [A = t\, a?\, [B = T2, 172], [D = T2, <j{\ and [E = t\,<j{\. 



to the conservation of world sheet current of space time energy-momentum of the test 
string, the total flux calculated around the path C must be zero, 



J A. 



(P;da-P°dT) = 0. (4.10) 

'ABDEA 

Note that one end of the string is attached to the boundary and the other end close to 
horizon is free. In the static gauge, Eq. (4.10) reduces to 

rti 



(4.11) 



t2 
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.^=0.80 
./i=1.20 



Figure 4. This plot shows the drag force from the analytic black hole solution as a function of T 
for chemical potential /i = 0.01, 0.80, and 1.20 respectively. Here we take v = 0.1. 




(a) (b) (c) 

Figure 5. This figure shows the drag force as a function of T for the chemical potential fi = 0.10 
(plot a), 0.34 (plot b), and 0.80 (plot c), respectively, in the hQCD model. Here the dashed curves 
stand for the behavior of drag force in confined phase which is denoted by CF in the figure. In the 
confined phase, in fact the drag force is not well defined, meaning that the dashed curves do not 
make any sense here. 



where p Xl is the x\ component of the total momentum at time t. Consequently the drag 
force is defined as 



dp Xl _ _ 1 £ 2 



2A a 



F *«» = It = = -^"ir"- (4 - 12) 

Finally we have to replace all the gravity parameters in terms of gauge theory parameters. 
Before doing that, we analyze the form of the constraints case by case. The exact forms of 
constraint for the solutions (3.9) and (3.15) are given respectively by 

v 2 = 1 - 4^1(3 sinh 4 (f^) + 2sinh 6 (f^)) + l^sinh 4 (f* 

Io V gW ( Io h 9(r)dr g(y) 5 rfy) dx ^ J£v g ^ dx 



It is always desirable to express the drag force in closed analytic form as a function of 
gauge theoretical variables. However it is very difficult to obtain analytic forms for the 
constraint (4.13) and the temperature (3.6). Instead we here solve them numerically and 
plot the drag force with respect to gauge theory parameters, e.g, temperature and chemical 
potential, so that the qualitative features of the drag force can be revealed. 

Certainly the analytic black hole solution (3.9) does not dual to a QCD model. As 
a warm-up exercise, we plot the drag force in Figure [4] for a model dual to this black 
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Figure 6. The figure shows the drag force as a function of T for three chemical potential fj, = 
0.10,0.34, and 0.80, respectively in the hQCD model. The dashed parts of the curves stand for the 
drag force in the confined phase which denoted by CF in the figure. Here we take v = 0.1. 

hole solution. We see that the drag force monotonically decreases with temperature and 
for a fixed temperature it becomes large as chemical potential grows. These features are 
qualitatively expected in realistic QCD. However, the behaviors of jet quenching parameter 
and screening length in the solution (3.9) are far away from QCD phenomenon and therefore 
we do not consider this solution from now on. 

In Figure [5] we plot the drag force for our hQCD model given by the solution (3.15) 
with different chemical potentials /x = 0.10, 0.34, and 0.80, respectively. We can see from 
the figure that for fixed chemical potential and temperature, the drag force increases with 
the velocity of the quark, while for fixed chemical potential and velocity, the drag force 
increases with temperature. These are expected features in QCD. In particular, let us 
note that in the low temperature region with small chemical potential, the drag force is a 
multi- valued function of temperature [see plot (a) and (b)], while it becomes a monotonic 
function with large chemical potential [see plot (c)]. This feature is closely related to the 
confinement /deconfinement phase transition in this hQCD model [15]. The dashed parts 
of curves in plot (a) and (b) denote the drag force in the confined phase and actually they 
do not make any sense here since drag force is not well-defined in the confined phase. Our 
result for the drag force in the deconfined phase is in agreement with the one in [54]. For 
comparison, in Figure [6] we plot the drag force versus temperature with three different 
chemical potentials \x = 0.10, 0.34 and 0.80, respectively. In this figure the velocity of 
quark is taken as v = 0.1. 

5 Jet Quenching parameter 

In this section, we use the AdS / CFT duality to compute the jet quenching parameter in our 
hQCD model. To employ the holographic principle in the computation of this parameter, 
the hQCD dual solution (3.15) as the bulk theory of gravity serves the purpose. The jet 
quenching parameter is related to the expectation value of light-like Wilson loop computed 
in adjoint representation. The gauge/gravity duality prescribes how to map the expectation 



-13- 



value of light-like Wilson loop in fundamental representation (W F (C)} into the exponential 
of the regularized extremal surface in the bulk with a boundary contour C located at z = 0, 
i.e., the AdS boundary, 

(W F (C))=eMS(C)]. (5.1) 

In the planar limit, considering the fact TrAdj = Tr 2 Fund , the relation between W F (C) and 
W A (C) can be easily established as 

(W A (C)) = (W F (C)f. (5.2) 

Now we start with the background black hole solution in string frame 

£ 2 e 2As ( dz 2 \ 

ds 2 s = ^ ( -f(z)dt 2 + — — + dx 1 dx 1 + dx 2 dx 2 + dx 3 dx 3 \ . (5.3) 

By introducing the light cone coordinates defined as 

the black hole metric (5.3) can be rewritten as 

p2 2A a /(i _ f(~\) „ „ 

ds 2 s = ^ J±2i((dx + f + (dx-f) - (1 + f(z))(dx + dx-) 

dz 2 

+J7-T + dx 2 dx 2 + dx 3 dx 3 ). (5.5) 

We take the gauge with r = x~(0 < x~ < L~), a = x 2 (—^- < X2 < -rr), and set the pair 
of quarks at x 2 = iy on x + = constant, x 3 = constant plane. In the limit with L~ S> L<i 
the string profile is completely specified by z = z(a). Following [33], using (5.1), (5.2), one 
has 

(W A (C)} = eM-^qL-L 2 ), (5.6) 



where the jet quenching parameter is defined as 

A 8y/2(S - S ) 
q = L-i 2 



(5.7) 



where S is the Nambu-Goto action of the string and So is the self energy from the mass of 
two quarks. 

Substituting the induced metric of the fundamental string into the Nambu-Goto action 
(4.1), we get 



2^7 / dTdayS-detg^ 



02 „2A S I '2 
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Since the integrand in (5.8) does not explicitly depend on a, one can regard a as time and 
the integrand as a Lagrangian. In this case the corresponding Hamiltonian is conserved. 
That is, we can have 

7r z'-£ = E, (5.9) 
dz 

where E is a constant and C is the integrand in (5.8). From this relation we obtain the 
equation of motion for z as 

2 ^ /W( ^™_ 1) . (5 ,0) 

We choose the boundary conditions z(±-y-) = and z (0) = 0. In that case, the 
turning point zt is determined by solving Eq.(5.10). Since z (a) is a real function, so the 
square of it should be non-negative. The realization of boundary condition z (0) = at 
the turning point requires the proper choices of zeros and the positivity region of the right 
hand side of Eq.(5.10). From the boundary conditions of the black hole solution 

lim/(z) = 0, lim/(z) = l, (5.11) 

z^Zh z— >0 

together with the fact that we are interested in the case with small E, it is clear that the 
factor " 4As (i-w)-e 2 * 4 i s always positive near the black hole horizon and negative near the 
boundary. To remove the region with a negative z 2 , we consider a modified boundary at 
z = 5. We assume that at z = z m i n , 



[ /lmm ^-l = 0, (5.12) 



z, 



L„ E 2 



nan 



and 5 > z m i n . In the region 5 < z < z^, thus, the factor {^S- — ~ 1] is always positive. 
So only viable solution of z 2 = is 

f(z) = 0^z T = z h . (5.13) 

That is, the turning point is just at the horizon. The distance between two quarks can be 
determined by 

L 2 r h . E 



— = / dz , — (5.14) 

As we are interested in the small L 2 limit, considering the smallness of E, we can expand 
Eq.(5.14) in terms of E as 

Lo f Zh z 2 e~ 2As E 2 f Zh e~ 6As z 6 

±1= dz z e + fL / dz e z +Q(E 4 ). (5.15) 

Inverting (5.15) suitably, we can obtain E up to the leading order of L 2 as 



- 15 - 



Thus we can obtain the string action 



L~£ 2 r Zh e 4As (l - f) 

1 dz [ J) . (5.17) 



s/2-kcI Js z 2 ^f{e AA °{l-f)-z A E 2 ) 

Clearly this action is divergent. The divergence comes from the contribution of mass of 
two quarks. With the gauge x~ = r and z = a, the self energy of two free quarks reads 




S = ^^ fa J^—LL m (5.1 8 ) 



Thus the regularized action up to the leading order of L2 is given by 

L-L\e 

&y/2na' //" dz 



= L ~3 e , 7Z. ■ \^a s + 0{L^). (5.19) 



With the definition of the jet quenching parameter (5.7), we finally reach 



(5.20) 



In fact the cutoff here can be removed by noting the fact that the integrand is regular 
inside the region < z < Zf L , i.e, from the horizon to the real boundary, 

z h z 2 e ~ 2A = r z h z 2 e~ 2As f 5 z 2 e~ 2As 



dz ; = / dz - - / dz ; (5.21) 

The second integral in the right hand side of the above equation smoothly vanishes in the 
limit (5—7-0. So the final expression for the jet-quenching parameter is 

e i . . 

=nr- (5-22) 



ira 



Because the black hole metric is still too complicated to obtain an analytical expression of 
the jet-quenching parameter in terms of physical parameters, we plot in Figure [7] the jet- 
quenching parameter as a function of temperature in the hQCD model with three chemical 
potentials \i = 0.10, 0.34 and 0.80, respectively. For large /i> p, c cases, the jet-quenching 
parameter decreases monotonically with temperature, which agrees with the one in [56] 
qualitatively. On the other hand, when fj, < fx c , the jet-quenching parameter is a multi- 
valued function of temperature in low temperature region and it decreases monotonically 
with respect to temperature in high temperature region. The multi-valued behavior of the 
jet-quenching parameter in low temperature region is clearly related to the first order phase 
transition between hadron phase (confined phase) and QGP phase (deconfined phase). The 
jet-quenching parameter confirms the hydrodynamical description of QGP phase and agrees 
with the real QCD expectation in high temperature. Once again, as the drag force in the 
confined phase, the dashed parts of curves in Figure [7] denote the jet-quenching parameter 
in the confined phase and thus they do not make any sense. 
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Figure 7. The figure shows the jet-quenching parameter as a function of T for three chemical 
potentials /i = 0.10,0.34, and 0.8, respectively, in the hQCD model. The dashed parts of curves 
stand for the jet-quenching parameter in the confined phase which are denoted by CF in the figure. 

6 Hot plasma wind and screening length 

The screening length is defined as the maximum length achieved by a quark-antiquark 
bound state at temperature T > T c , beyond which the pair dissociates. For quark-antiquark 
pair, the energetically favorable configuration in the dual gravity theory is a fundamental 
string with both ends attached to the boundary. The attached endpoints correspond to 
the qq pair whereas being separated beyond the screening length, thus dissociated from 
each other, the pair maps into two separate strings hanging from the boundary. In [55], 
the screening length is computed in the rest frame of qq pair and the plasma wind flows 
at a constant speed v for the hot J\f = 4 SYM plasma. This setup is identified with a 
quark-antiquark pair moving in hot N = 4 SYM plasma. In this section, we compute the 
screening length for the hQCD model (3.15) in the same way as in [55]. 

In the static frame of qq pair, we assume that the hot plasma is moving with velocity 
v in the negative X3 direction. The Wilson loop we are interested in lies in the t — x\ plane 
specified by the length T and L respectively. We assume T ^> L such that the string world 
sheet is invariant under translation along the time direction. The boost we are considering 
is defined as 




where cosh 77 = 7, sinh r] = jv and 7 = 1/^1 — v 2 is the Lorentz boost factor. With the 
Lorentz transformation, we obtain the boosted black hole metric in string frame 
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with the following boundary conditions 

z(a = ±—) = 0, z(a = 0) = z c , z {a = 0) = 0. (6.4) 
Thus the world sheet metric induced on the boosted background is given as 

g TT = - J ff(z)(l-(l-/)cosh 2 r ? ), 
9tct 9<jt 0) 

g aa = H(z)[l + (l + j---)]. (6.5) 
Then the Nambu-Goto action for the string takes the form as 



S = ~—, 1 2 daH(z)\ (l-(l-f) cosh 2 r ? )(l + ^-). (6.6) 

As the Lagrangian C in (6.6) does not depend on a explicitly, the corresponding Hamilto- 
nian is conserved and can be viewed as a constant of motion 

= (6-7) 
With this we can cast the equation of motion in the form as 



//[#2(l-(l-/)cosh^)-g2] / 
z' = . (6.8) 

It is evident from the constraint (6.8) that at the horizon, z = Zh, where f(zh) = 0, the 
factor ^-(1 — (1 — /) cosh 2 rj) — 1 = — ^ sinh 2 77 — 1 is always negative. At the boundary, 

/(0) = 1, the factor ^(1 — (1 — /) cosh 2 rj) — 1 = ^ — 1 is always positive for small values 
of q < H. Therefore in the range < z < z^ there must be a location {z = z c ) where 
— (1 — /) cosh 2 rf) — 1 switches its sign. Accordingly z = z c is the physical turning 
point of the string configuration. The string can not be stretched up to the horizon as z' 
is an imaginary quantity in the region z c < z < Zh- By solving the equation 

f(z c )H 2 {z c ) cosh 2 7] _ H 2 {z c ) sinh 2 V _ 1 = Q ^ Q , 

q 2 q 2 

the turning point can be numerically determined. Then one can obtain the binding energy 
between the quark and antiquark pair through calculating the action (6.6) with constraint 
(6.8) 

V = -^, (6.10) 



where Sq is given by 



T 



So = r / dzyJ-G tt G zz . (6.11) 
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Figure 8. Plot (a) shows the quark- ant iquark distance as a function of q for a fixed rapidity, while 
plot (b) shows the binding energy with respect to the distance. In both plots we fix the chemical 
potential yU = 0.10,0.34, and 0.80, respectively. We have set a same temperature T to obtain these 
curves. 



The distance between quark and antiquark can be calculated from (6.8) as 

1 

dz . (6.12) 

H^f f[(l - cosh 2 V (l-f))-£] 

It is not possible to work out the integration in (6.12) explicitly. To determine the screening 
length, we plot the distance L with respect to the constant of motion q for a fixed rapidity 
7] in Figure [8] (see plot (a)). It turns out that for a fixed value of rapidity, there exists 
a maximum for L, which is regarded as the screening length L s = L max (r]) / (irT) . Plot 
(b) in Figure [8] shows the binding energy V given by (6.10) with respect to L. One can 
see from plot (a) that the quark-antiquark distance starts from zero when q is also zero, it 
increases sharply with respect to q, reaches its maximum at a certain q, and then decreases 
monotonically to zero at some finite q. In between these two zeros, there exists a single 
L = L max beyond which there is no solution of Eq.(6.12). This implies the quark-antiquark 
pair dissociates beyond L = L max . We identify L max (r}) / (nT) with the screening length 
L s . For the fi = 0.1 case, L max ~ 1.4 and L s ~ 1.4/(-7rT) ~ 0.45/T, close to the lattice 
calculation L s ~ 0.5/T [57] of the static potential between heavy quark and antiquark in 
QCD. Plot (b) shows that there are two branches for the binding energy in the region 
L < L max . The branch with dashed curves has a higher energy than the one with solid 
curves. This implies that the branch with dashed curves is physically disfavored. 

The screening length L s (rf) as a function of rapidity can be obtained numerically as 
illustrated in Figure [9] . One finds that it decreases with velocity which indicates that the 
quark-antiquark pair dissociates at a lower temperature as it is moving. This behavior is 
also observed in [54]. If the qualitative behavior holds for QCD, it will have the consequence 
for quarkonium suppression in heavy ion collision. Additionally, our results show that the 
case with smaller chemical potential has a larger screening length. 




- 19 - 



Z™, CoshM 1 ' 




Figure 9. The screening length versus the rapidity rj for the cases u = 0.1,0.34, and 0.8, respec- 
tively. We have set a same temperature T to obtain these curves. 



7 Free energy and phase transition 

In this section, we would like to calculate the free energy of the black hole solution dual to 
the hQCD model by following [46]. Ref. [46] studies the free energy of a generic Einstein- 
Maxwell-Dilaton system. The system considered in [46] is given by 



-9 



R 



^F 2 -Ud<$>) 2 + V{<$>) 



+ Igh 



I GH = 2Mfr 1 J dM d d xVhK 
with following ansatz 



(7.1) 



da 2 



-f{u)dt z + dx'dx 1 + 



du 2 
/(«) 



, A = A t (u)dt , $ = 



(7.2) 



where Mp is the Planck mass, and K is the extrinsic curvature of the finite boundary dM 
with induced metric h. By the following transformations 



A. = e -A{u{z))j_ = e -2A(u(z)) ( _f_ _ dA{u{z)) d 

du dz' du 2 \dz 2 dz dz 



f(u(z)) = f(z), 



$(u) 



Z(fb) = 1, d = 4, Ml" 1 



1 



I(u(z)) = A e (z) - log(z), - V{u(z)) = V E (z), 



(7.3) 



we can change the system and ansatz to ours discussed in this paper. Here we have chosen 
A(u(z)) = A e (z) — iog(z) as a gauge. Varying the action (7.1) yields equations of motion 
for the gravitational field 



y F w FP-jg^F pa F^ =0(7.4) 
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The Maxwell equation and the explicit forms of Einstein equations read 

I ( e{d - 2,AzA <) ~ ■ (7 ' 5a) 

2(d- 1)1 + 6 2 = , (7.5b) 
/ + dAf - e~ 2A ZA 2 = , (7.5c) 

(d - 1)1/ + (d(d -l)A - 1 -&\f-V+ \ze- 2A A 2 t = , (7.5d) 

where an over dot stands for the derivative with respect to u. Eq.(7.5a) is the Maxwell 
equation and the other three Eq.(7.5b) (7.5c) (7.5d) are obtained from Einstein equations. 
One can easily check that Eq.(7.5b) and Eq.(7.5c) correspond to Eq.(2.8) and Eq.(2.9), 
respectively. In addition, the equation of motion for the scalar field is 



--d u (V=gd u <t>) + V'($) = °. ( 7 - 6 ) 



where the prime denotes the derivative with respect to $. From (7.5a), one obtains 

A t (u) = u+I du, (7.7) 

where fj, and p are the chemical potential and charge density of the black hole solution, 
respectively, and no stands for the UV boundary. 

By defining the superpotential W in the following way [46] 

$ = W'($), (7.8) 

the equation (7.5b) can be solved as 

A = -^L (79) 
2(d-l)' 1 ' 

Equivalently, 

1 /"* ~ W($) /"* 8$ 

A($) = A --—— 6$—y-, u = u + — -j,-, (7.10) 



2{d - 1) J$ W'$) ' U W'{&) 

where A = A(<f> ), $(u = uq) = <J>o- 

The temperature T associated with the black hole in metric (7.2) is given by 



47T 



e A f 



(7.11) 

U=Uh 



where Uh denotes the black hole horizon. The entropy density for the black hole solution 
(7.2) is 



S = ^e^ 1 ^ = 47re {d - 1)io -^*o h5l> ^, ($ h := $«)). (7.12) 
where M£~ l = 1 has been taken. 
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Figure 10. The free energy of the background black hole solution with respect to temperature. It 
shows that in the small p < p c case, the free energy has a multi- valued behavior, while it is absent 
in the large p > p c case. The parameters £ = 1, 167rG5 = l,g g = 1, k = 0.3GeV are taken here. 

The on-shell action / has a bulk contribution Ie and a boundary contribution Iqh 
from the Gibbons-Hawking term. Here we evaluate it with an UV cutoff at u = uq or 
z = e. A standard computation gives us with the action density 



With the equations (7.3) (7.5c) and A t (uh) = 0, the free energy density is found to be 



case d = 4 and the expression of the free energy is written in the ansatz (2.6) in Einstein 
frame. The first term in Eq.(7.14) corresponds to the black hole mass and the last term 
can be expressed as —\xp. 

With the analytical result (7.14), one can obtain the free energy of the black hole 
solution (3.15) with respect to temperature with various chemical potentials numerically. 
Figure [10] shows the free energy density of the black hole solution. It can be seen that 
the free energy density is a multi-valued function of temperature in the case p < p c and 
it becomes monotonically in the large p region p > p c . Clearly it indicates that in the 
small p region there exists a first order phase transition, while it is absent in the large p 
region. The existence of the first order phase transition in the small p region is consistent 
with the analysis of Wilson loop in this black hole background [15]. The existence of the 
critical point at p = p c is in agreement with recent lattice calculation given in [58] . Let us 
notice that in [15] we claim that there is a crossover or higher order phase transition in the 
large p > fi c region, while the free energy calculation here shows no such phase transition. 
These two conclusions are in fact not in contradiction with each other. The reason is as 
follows. The free energy calculation here only concerns with the black hole background, 
from the point of view of QCD, the black hole background does not include the degrees of 




(7.13) 



W($o) - TS + p 2 / ^ , 




Figure 11. The phase diagram of the hQCD model. The red and dashed red curves are obtained 
in [15] through studying VEV of Polyakov loop. The two curves stand for first order and continuous 
phase (or second order) transition respectively. The thick black curve is obtained by calculating 
the free energy in this paper. The black dot denotes the critical point [15]. 



freedom of quarks, while the Wilson loop calculation in [15] is related to the dynamics of 
quarks. 

In Figure [11] we plot the phase diagram of the hQCD model. The difference between 
the confinement / deconfiement phase transition lines come from the different consideration 
as mentioned above. 

At the end of this section, let us stress that the so-called Hawking-Page phase transition 
means the transition between the black hole in AdS space and thermal gas in AdS space [47] , 
while in the above analysis, our phase transition happens between small black hole and big 
black hole in the region with small chemical potential (M < fj, c , such a transition does not 
occur as fj, > \i c . It is therefore of some interest to investigate whether there exists any 
transition between black hole solution and thermal gas solution in our hQCD model. For 
this aim, we have to first find the thermal gas solution. Unfortunately, in this potential 
reconstruct approach, to find out the thermal gas solution is not an easy job since we have 
not an explicitly analytic form of the potential for the dilaton field in our hQCD model. But 
even so, we still can make some arguments which lead to the conclusion that there does not 
exist any phase transition between black hole and thermal gas solutions in the hQCD model. 
The main reason for this comes from the fact that in the hQCD model, for any temperature, 
there exist corresponding black hole solutions. This can be seen from Figure [2]. It shows 
in plot (a) that when fj, > fj, c , the black hole temperature increases monotonically from zero 
with the smallest black hole radius, which corresponds to an extremal black hole, while 
when fi < fie, plot (b) shows that as T < T(z m ), there exists a region with stable black 
holes with small horizon radius, including the extremal black hole; and as T > T(z m ), there 
exists the region with stable black holes with larger horizon radius. In addition, for a given 
temperature, it is generally believed that the black hole solution has much large entropy 
than the thermal gas, so that the black hole phase is dominated. As a result, whatever 
the chemical potential is, the black hole solutions are always dominated over the thermal 
gas solutions. This is quite different from the case of Schwarzschild-AdS black hole. In 
the latter case, there exists a minimal temperature, below which there is no black hole 
solution. This leads to the conclusion that in low temperature the thermal gas solution 
is dominated, while in high temperature the black hole solution is dominated, and thus 
there must exist the Hawking-Page phase transition between the black hole solution and 
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thermal gas solution at a certain temperature. Therefore in our case, black hole solutions 
are always dominated and there does not exist any phase transition between black hole 
solution and thermal gas solution. But as we analyzed above, there exists indeed the phase 
transition between small black hole and big black hole in our model with \x < \i c . This 
phase transition can also be understood as the confinement /deconfinement phase transition 
in QCD since clearly the small black hole has a much less entropy than the large one. 

8 Conclusion and discussion 

In this paper, we have continued to study a holographic QCD model proposed in [15], 
in an Einstein-Maxwell-Dilaton system. At first we have generalized case with a non- 
minimal coupling between Maxwell field and dilaton field, and given a generic formulism 
for generating a set of exact and asymptotic AdS black hole solutions in the EMD system. 
After briefly reviewing the main features of the hQCD model, we have studied some aspects 
of QGP phase of the hQCD model by calculating some quantities such as drag force, jet 
quenching parameter and screening length. The calculations show that the behaviors of 
those quantities are consistent with the expectation from real QCD. 

It is found that the drag force increases monotonically with temperature which is quite 
good consistent with real QCD phenomenon in the larger chemical potential region with 
fJ> > A*c- in the small chemical region with fj, < fi c , the drag force also monotonically 
increases in the high temperature region, while in the low temperature, it shows a multi- 
valued behavior. Note that in the case [i < /i c , the solution is dual to the confined phase 
of QCD. In that case, the drag force is not well defined. Therefore the change from 
the multi-valued behavior to the monotonic behavior just manifests the existence of the 
first order phase transition. The jet quenching parameter has monotonically decreasing 
behavior versus temperature, which is also consistent with QCD experiments in /i > fj, c . 
For the [i < fj, c case, the jet quenching parameter agrees with real QCD expectation in high 
temperature and once again, it shows the multi-valued behavior in the low-temperature 
region. As in the case of drag force, the multi- valued behavior of jet quenching parameter in 
the low temperature region is consistent with the existence of first order phase transition in 
this hQCD model. For the screening length we have plotted the separation between quark 
and anti-quark with respect to the constant of motion q. It is clear from the plot that for 
both cases when [i > [i c and /U < \x c the dipole dissociates beyond a maximum separation 
distance, namely the screening length L s . We have also calculated the binding energy as 
a function of separation distance. In addition, we have presented L s (rj) and found that 
there are qualitative consequences for quarkonium suppression in heavy ion collisions in 
this hQCD model. 

Finally, we have calculated the free energy of the background black hole solution and 
further confirmed that there exists a first order phase transition in small region between 
small black hole and big black hole [15]. When /i > /i c , the phase transition is absent from 
the point of view of the free energy. The existence of the critical point is consistent with 
the result in [15]. Further we have argued that there does not exist any phase transition 
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between black hole solutions and thermal gas solutions in this model, and the main reason 
is given. 

In this work we have studied some aspects of QGP phase in the hQCD model proposed 
in [15]. The results are encouraged and are consistent with the expectation of real QCD. 
Thus it would be of great interest to further investigate the hQCD model. For example, 
it is required in the model to study the spectra of hadrons, chiral phase transition and 
its phase diagram [59] [60] [61], hydrodynamical properties of QGP, and color flavor locked 
phase [62] [63] [64], etc. Furthermore, quantum corrections [65] to above physical quantities 
are also deserved to consider in the coming works. 
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